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We propose a hybrid quantum architecture for engineering a photonic Mott insulator-superfiuid phase transi¬ 
tion in a 2D square lattice of superconducting transmission line resonator (TLR) coupled to a single nitrogen- 
vacancy (NV) center encircled by a persistent current qubit. The localization-delocalization transition results 
from the interplay between the on-site repulsion and the nonlocal tunneling. The phase boundary in the case 
of photon hopping with real-valued and complex-valued amplitudes can be obtained using the mean-field ap¬ 
proach. Also, the quantum jump technique is employed to describe the phase diagram when the dissipative 
effects are considered. The unique feature of our architecture is the good tunability of elfective on-site repul¬ 
sion and photon-hopping rate, and the local statistical property of TLRs which can be analyzed readily using 
present microwave techniques. Our work opens new perspectives in quantum simulation of condensed-matter 
and many-body physics using a hybrid spin circuit QED system. The experimental challenges are realizable 
using current available technologies. 

PACS numbers: 03.67.Lx, 05.30.Rt, 42.50.Ct 


I. INTRODUCTION 


The microscopic properties of strongly correlated many- 
particle systems emerging in solid-state physics are in general 
very hard to access experimentalljnil. So how to simulate the 
properties of condensed-matter models using nontraditional 
controllable systems is highly desirable. Recently, the investi¬ 
gation of quantum simulation in the photon-based many-body 
physics has received much attention in different systems^Hl 
Especially, there has been a great interest in mimicking quan¬ 
tum phase transition (QPT) of light with scalable coupled res¬ 
onator array in the context of cavity/circuit quantum electro¬ 
dynamics (QEDjP®, which provides a convenient controllable 
platform for studying the strongly correlated states of light 
via photonic processes. On the other hand, the artificially en¬ 
gineered hybrid devices ca n perm it measurement access with 
unique experimental control^®!, and it is intriguing to employ 
a well-controllable quantum system with a tunable Hamilto¬ 
nian to simulate the physics of another system of interest. This 
paradigm has promoted many experimental/theoretical pro¬ 
posals on probing the light phase and opened various possi¬ 
bilities for the simulation of many-body physics. 

In this work, we develop an optical system for engineering 
the strongly correlated effects of light in a hybrid solid-state 
system. We consider a 2D square lattice o f coup led tlrP, 
where each TLR is coupled to a single Nencircled by 
a persistent current qubit (PCQ). We show that the compe¬ 
tition between the NV-PCQ-TLR interaction and the nonlo¬ 
cal hopping induces the photonic localization-delocalization 
transition. Subsequently the Mott insulator (MI) phase and 
the superfluid (SE) phase can appear in a controllable way. 
The phase boundary in the case of photon hopping with 
real/complex-valued amplitudes can be obtained using the 
mean-field approach. Also, the quantum jump technique is 
employed to describe the phase diagram when the dissipation 
is considered. Einally, the possibility of observation of the 
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FIG. 1. (Color online) (a) Schematic circuit for the resonator lattice, 
where each TLR is coupled to a single NV center encircled by a 
PCQ, and the circles denote the central coupler, (b) The subsystem 
consisting of NV and PCQ, where NV is at the center of the loop. 
The PCQ is made up of three Josephson junctions, and it couples to 
the NV via the magnetic field at the center of the loop generated by 
the persistent currents in the loop. The energy diagram of the NV is 
shown in the red box. 


QPT is discussed by employing experimentally accessible pa¬ 
rameters. 

In our architecture, one can tune independently the on-site 
emitter-field interaction and the nonlocal photonic hopping 
between adjacent TLRs. This permits us to systematically 
study the localization-delocalization transition of light in a 
complete parameter space. The main motivation for building 
such a hybrid system is to combine several advantages: in situ 
tunability of circuit elementJ^ spectroscopic technology for 
state readout, peculiar characteristics of NV (e.g., individual 
addressing and long cohere nce t ime at room-temperatur^^^, 
and scalability of TLR arrayP^^tSl. Recently, D. Underwood 
et al experimentally fabricated 25 arrays of TLRs and demon¬ 
strated the feasibility of quantum simulation in circuit QED 
systerrP^. E. Lucero et al experimentally characterized a 
complex circuit composed of four phas^ubits and five TLRs 
to realize intricate quantum algorithms^. The progress ren- 
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FIG. 2. (Color online) The dependence of effective on-site repulsion 
U± on the photon number n under the different detunings A, where 
the solid (dotted) line denotes U+ ((/_). 

ders the TLR lattice a good platform for studying condensed- 
matter physics with photons and makes our scheme to be more 
practical. 

II. MODEL 

As illustrated in Fig. 1, we consider a 2D lattice of coupled 
TLRs, where the basic unit consists of a TLR coupled to a 
single NV encircled by a PCQ, which acts as an interconnect 
to greatly magnify the NV-TLR coupling by several orders of 
magnitude, compared with the direct NV-TLR coupling (far 
below the linewidth of TLR with dozens o f kHz ) resulting 
from the vacuum fluctuations of the photonJ^^^. The TLR 
is made of a superconductor line interrupted by two capaci¬ 
tors at its ends. In the microwave domain, it can be treated 
as a quantum LC harmonic oscillator, = C 0 r{a^a + 1 /2) 
{Ji = 1), where w,. = yjljLrCr is the corresponding eigen- 
frequency with inductance Lr and capacitance Cr- The PCQ 
located at an antinode of TLR’s magnetic field is formed 
by a superconducting loop interrupted by three Josephson 
junctionJ^. When the loop is biased by half a magnetic flux 
quantum, the device is an effective two-level qubit made up 
of two countercirculating persistent currents with the Hamil¬ 
tonian - The NV can be modeled as a three- 

level system in the triplet ground-state subspace consisting of 
pA,mj = 0^ and = +1^ The Hamiltonian is = 

+ D{S^ - 2/3), where is the electronic gyromag- 
netic ratio, DI2n ~ IXlGHz is the zero-field splitting, is a 
perpendicular magnetic field at the NV and 5^ is the spin-l-z 
operator. 

The PCQ magnetically couples to TLR via mutual induc¬ 
tance, - -fi ■ B, where fi is the magnetic dipole of 

PCQ induced by the persistent circulating currents and B is 
the magnetic field at PCQ induced by the current in the cen¬ 
tral conductor of TLR. When cj^ ~ ojq, we have - 

gia^cr^ + acr^) after rotating wave approximation, where g - 
{Ipfxor^ld) Vwr/2Lr, r (Ip) is the radius (persistent circulating 
current) of the PCQ loop, and d is the distance between PCQ 
and central conductor of TLR. The sizable changes of mag¬ 
netic flux within the loop induced by Ip presented in the PCQ 


lead t o sm all shifts in the transition frequencies (nis —> +1) 
of n\ 122| 23 [ Through this small change in magnetic field the 
PCQ can couple to the NV via Zeeman term, — ^tjct^S 
where rj = IpHoyJr. 

The basic unit in our system is thus governed by the Hamil¬ 
tonian Hi = H- + fj^v + fjT-P + fjN-P 

photonic tunneling in our model can be realized by a central 
couple]!^ which serves as individual tunable quantum trans¬ 
ducers to transfer photonic states between adjacent TLRs. We 
have presented a new paradigm for 2D TLR lattice coupled 
to solid-state spins. We have shown that specially engineered 
resonator lattice provides a practical platform to couple both 
individual spin and superconducting qubit, and engineer their 
interactions in a way that surpasses the limitations of current 
technologies. This can provide new insights to many-body 
physics. 

III. MOTT-SUPERFLUID TRANSITION 

We study the full Hamiltonian of the 2D square lattice by 
adding the on-site chemical potential and the nonlocal mi¬ 
crowave photon hopping between adjacent sites. The Hamil¬ 
tonian is given by 

H — ^ Up + ^ k(pYO.paq — ^ lipNp, ( 1 ) 

p (p,q) P 

where k^p^q) - 2ZoC(p,^)(Wr + dp)(cJr + 5q) are photonic tun¬ 
neling rates between resonators p and q, which are set by the 
tunable mutual capacitance between resonator ends with 
characteristic impedance Zq and frequency shift dp due to ran¬ 
dom disorders Since Ur » 6p, one can assume that k^p^ - 
k - 2ZoC<jj^ without disorder for nearest-neighbor resonators, 
and k(^p,q) - 0 for other resonator pairs, pp is the chemical po¬ 
tential at the p-th site. The conserved quantity in our system 
is the total number of excitations Np - a^Op + crlcTp + p 
with S' (cr') the spin-1 (-1/2) operators (i - x,y, +). 

The photon-dependent eigenstates of the Hamiltonian 
is dressed states \ms, +, n) with the eigenvalues E\m„±,n) - 
na)r + (-A + 2nA± yjAng^ -H + 2Aq + q^)l2 +x(fns), where 
Xints) - D(3 X 1'"' -2) + niYeBz is the eigenenergy of 
Here A = cUr - wq is the detuning, n is the number of ex¬ 
citations in the resonator and |±) = (\e) ± \g))/ V2. In our 
case, the dynamics is governed by the Jaynes-Cummings (JC) 
type of interaction, which enables the interconversion of qubit 
excitations and photons, and provides the effective on-site re¬ 
pulsion. Meanwhile, pairs of TLRs are coupled by the two- 
site Hubbard model via one-photon hopping. The difference 
between the Bose-Hubbard model (BHM)^ and our model 
is that the conserved particles are the polaritons rather than 
the pure bosons in BHM, and the effective on-site repulsion 
U±(n) - £’|m,,±,«+i> - E\m„±,n) decreases with the growth of 
photon number, and U±(ri) ^ 0 in the limit of large n and 
A = 0, as shown in Fig. 2, while it is a constant in BHM. 

The phase diagrams can be distinguished using the cor¬ 
responding order parameters. Here we choose the SF or¬ 
der parameter ^ - {Up) (set to be real) to differentiate be¬ 
tween insulator-like and SF-like states. Using the mean-field 
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FIG. 3. (color online). The phase diagrams m ^ ~ k plane for dif¬ 
ferent sets of NV-PCQ coupling rj and tunable magnetic fields 
applied on NV, where ji and k are the chemical potential and the 
photon hopping rate. The common parameters are g = I, air = 200, 
D = 100 and = -10^. The other parameters are set as (a) 77 = 0.01, 
= O.OOOST and A = g, the phase boundaries are plotted in (d), 
where A = 2g (solid), g (dashed) and 0 (dot-dashed); (b) jj = 1.2, 
B^ = -0.3 and A = 0, the phase boundaries are plotted in (e), where 
B, = -0.3r (solid), O.OOOST (dashed) and 0.3r (dot-dashed); (c) 
T] = 0.75, Sj = 0.3r and A = 0, the phase boundaries are plotted 
in (f), where 77 = 0.01 (solid), 77 = 0.75 (dashed) and 77 = 1.5 (dot- 
dashed). 


theorjEH we decouple the hopping term as + 

Op - (up^ {^q)^ the resulting mean-held Hamiltonian can 
then be written as a sum over single sites, 

rMF ^ ^ [^0 _ ^ ^ ^ (2) 

p 

where z = 4 is the number of nearest neighbours. Noting that 
,S^ - 0 , therefore, we can treat in the mean-field 
Hamiltonian as a c-number and — +1,0. Minimizing the 
ground state energy of the Hamiltonian with respect to \j/ 

for different values of 72 and k, we obtain the mean field phase 
diagram/boundary in the ( 72 , k) plane when 77 varies from the 
weak coupling regime (77 <s: g) to the strong coupling regime 
(77 > g) under the resonant/detuning case. The features of 
Fig. 3 are rich. The regime where ijj - 0 corresponds to the 
stable and incompressible MI lobes characterized by a fixed 
number of excitations at per site with no variance. In each MI 
lobe, due to the nonlinearity and anharmonicity in the spec¬ 
trum originating from the photon blockade effecpni, the strong 
emitter-field interaction leads to an effective large polariton- 
polariton repulsion which freezes out hopping and localizes 
polaritons at individual lattice sites. By contrast, strong hop- 



FIG. 4. (color online), (a) The order parameter in 77 ~ k plane. The 
parameters are g = 1, ut,. = 200, D = 100, = -10'^, 77 = 0.01, B^ = 

0.0005r, a = 0.2, and A = 0. The corresponding phase boundaries 
are plotted in (b), where the solid, dashed, and dot-dashed line denote 
a = 0, 0.2, and 0.24, respectively. 


ping favours delocalization and condensation of the particles 
into the zero-momentum state, namely, i/r 4 0 indicates a SF 
compressible phase with the stable ground state at each site 
corresponding to a coherent state of excitations. 

Analogous to the BHM, the physical picture behind is 
that the MI-SF phase transition results from the interplay be¬ 
tween polariton delocalization and on-site repulsive interac¬ 
tion. Therefore, the phase boundary primarily depends on the 
ratio of the photon-hopping rate to the on-site repulsion rate. 
When the on-site repulsion dominates over hopping, the sys¬ 
tem should be in the MI phase, otherwise the system will be 
in the SF phase. From the expression of the parameter 77 and 
g, we can find that reduction of the size of the PCQ loop will 
increase 77 but decrease g, and the adjustment of the distance d 
only affects TLR-PCQ interaction. Furthermore, the detuning 
A is also tunable by varying the magnetic field applied on NV. 
In Fig. 3, one can find that the size of the MI lobes varies with 
A, with the largest Mott lobes found on resonance. 

Further insight to the transition can be gained when the pho¬ 
ton hopping with complex-valued amplitude exists in Eq. (1), 
where the hopping process becomes - Yj(p,q} ^(p,q)^“^^'’'‘'^‘^p^q 
with <p(p,q) = -<P(q,p) we set k(p^q) - k. We emphasize that 
this process is possible if the intermediat e coup ling elements 
are used to break time-reversal symmetrjEE^El. The param¬ 
eter k(p^q)e‘'^'-i’-‘’'> provides a new regime in the dynamical evo¬ 
lution of the system. The sum of the tunneling phases along 
a closed loop surrounding the plaquette is 2 ( 0 ^,+^^ + (f>p^q+i - 
4’p,q ~ 4’p+i,q+i) - 22 ra', which is actually the flux quanta per 
plaquette. Assuming that as are all equal, the total Hamilto¬ 
nian under mean-field approximation reads 

= '^p\-f^p - zkiff cos(2na)(ap + Op) 

+zk>j?' cos(27ra) - gipNp]. (3) 

The results are exhibited in Fig. 4, we find that the boundary 
line gradually shifts to the right as a enhances in the interval 
[0,1/4]. Because of the spatial variation of tunneling phase, 
the wave function of a polariton from one lattice site to another 
acquires a nontrivial phase (Aharonov-Bohm phasewhich 
actually reduces the effective hopping rates. 
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FIG. 5. (color online). The order parameter ij/ m ~ k plane under 
the different dissipative rates (a) F = k = 0.01; (b) F = k = 0.05; 
(c) F = K = 0.1; (d) F = K = 0.15. The other parameters are g = 1, 
ojr = 200, D = 100, je = -10\ r] = 0.01, B, = 0.00057’, and A = 0. 


IV. DISSIPATIVE EFFECTS 


Generally, nonequilibrium processes such as dissipative ef¬ 
fect, are crucial in solid-state devices. We show that the signa¬ 
ture of the localization-delocalization transition remains even 
in the presence of the engineered dissipation by the quantum 
trajectory method^. The non-Hermitian Hamiltonian is for¬ 
mulated as 


jj — jjMF 

tide - ti - — 


/r 


- - 


iK 

2 2-ii. 


apQp, 


(4) 


where k = 4ZyC^„,<u^ is the decay rate of TLR, and T is the 
decay rate from the effective excited state \e) of PCQ. In our 
case, the dissipative effects result from the unavoidable inter¬ 
action between the PCQ/TLR and the corresponding Marko¬ 
vian environment, for example, the interaction between the 
output of the TLR and the corresponding vacuum field will 
result in a photonic escape rate with k to the continuum. Here 
the dissipative effects of NV are negligible, compared with 
K and r. The phase diagrams under dissipative effects are 
displayed in Fig. 5. Once the hopping rate is increased be¬ 
yond a critical value, the system is expected to undergo a non¬ 
equilibrium QPT from a MI phase, where the initial photon 
population is self-trapped, to a SF phase with the dynamical 
photon population imbalance coherently oscillating between 
pairs of tlrP. Furthermore, another obvious feature is that 
the size of MI phase becomes larger as the growth of dissi¬ 
pative rates. Note that the effective nonlinearity at per site 
becomes stronger at lower exciton numbers, which implies 
that the dissipative effect (inducing the decrease of the exciton 
numbers) favours the MI phase. As a result, the dissipation re¬ 
sults in the dynamical switching from SF phase to MI phase 
and causes the increment of the size of MI phase. 


V. EXPERIMENTAL FEASIBILITY 

Firstly, we briefly stress the relevant experimental progress. 
From theoretical standpoint, it is possible to fabricate large 
arrays to observe many-body physics of interacting polaritons 
since resonators and qubits can be made lithographic alljEl. 
Actually, it is indeed feasible to couple over 200 (or 1000) 
TLRs with negligible disorders (on the order of a few parts 
in lO”*) in a 2D lattice using a 32mm x 32mm sample or a 
full two-inch wafei®. Secondly, how to probe quantum many- 
body states of light is still an open que stion in photonic quan¬ 
tum simulatiorPS. The previous workfl^ suggested to mea¬ 
sure the individual TLR through mapping the excitations into 
the qubit followed by state-selective resonance fluorescence 
spectrum, but a remaining technical challenge is the realiza¬ 
tion of high-efflciency photon detectors. Alternatively, the 
local statistical property of TLR can be analyzed readily us¬ 
ing com bined techniques of photon-number-dependent qubit 
transitiorPIEll and fast readout of the qubit state through a sep¬ 
arate low-Q resonator modd^, for which the high-efficiency 
photon detectors are not required. Experimentally, transmis¬ 
sion and reflection measurements for circuit QED atTays have 
been impleme nted successfully in small system with one or 
two resonator^212S) Therefore, in order to distinguish be¬ 
tween different phases of the system, one can also experi¬ 
mentally probe beyond transmission, such as two-tone spec¬ 
troscopy and second-order coherence function (photon statis¬ 
tics) to reveal additional information. The tunability of cou¬ 
pling strengths in our system enables one to measure these 
quantities relatively straightforwardly. 

Einally, we survey the relevant experimental parameters. 
Given the flexibility of circuit QED, we can access a wide 
range of tunable experimental parameters for TLR-PCQ cou¬ 
pling strength g and hopping rate k(p,q). Taking - 2 nH, 
(i)rl2n - 6 GHz, Ip - 600 nA, and r - 0.2 yum, we get 
ril2n ^ 140 KHz and gl2n ~ [1.8,180] MHz when the dis¬ 
tance d varies from 5 jim to 50 nm. Eurthermore, the hopping 
rate k(p,q) depends on the tunable mutual capacitance C(p^q) be¬ 
tween resonator ends. In Ref®! the authors measured devices 
with photon hopping rates kl2n form 0.8 MHz to 31 MHz in 
resonators lattices. On the other hand, the electron-spin relax¬ 
ation time T i of NV ranges from 6 ms at room temperaturd^ 
to 28 ~ 265 s at low temperaturd^. In addition, later experi¬ 
mental progresJ^with isotopically pure diamond has demon¬ 
strated a longer dephasing time to be r 2 = 1.8 ms. Therefore, 
the dissipation and decoherence of NV are negligible. 


VI. CONCLUSION 

We have devised a concrete hybrid system to engineer pho¬ 
tonic MI-SE phase transition in a 2D square lattice of TLRs 
coupled to a single NV encircled by a PCQ. We And that the 
interplay between the on-site repulsion and the nonlocal tun¬ 
neling leads to the photonic localization-delocalization tran¬ 
sition. In the presence of dissipation, the phase boundary can 
be obtained by the mean-field approach and the quantum jump 
technique. Eacilitated by high levels of connectivity in circuit 
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QED, experiments combining both scalability and long coher¬ 
ence times are expected in the coming few years, at that stage 
the investigation of photonic QPT using TLR lattice systems 
can therefore be easier to realize. 
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